The Dynamic Structure Factor of the ID Bose Gas near the Tonks-Girardeau Limit 
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While the ID Bose gas appears to exhibit superfluid response under certain conditions, it fails 
the Landau criterion according to the elementary excitation spectrum calculated by Lieb. The 
apparent riddle is solved by calculating the dynamic structure factor of the Lieb-Liniger ID Bose 
gas. A pseudopotential Hamiltonian in the fermionic representation is used to derive a Hartree-Fock 
operator, which turns out to be well-behaved and local. The Random-Phase approximation for the 
dynamic structure factor based on this derivation is calculated analytically and is expected to be 
valid at least up to first order in I/7, where 7 is the dimensionless interaction strength of the model. 
The dynamic structure factor in this approximation clearly indicates a crossover behavior from the 
non-superfluid Tonks to the superfluid weakly-interacting regime, which should be observable by 
Bragg scattering in current experiments. 
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The emergence of superfluidity at low temperatures 
is one of the most dramatic manifestations of quantum 
many-body physics in nature. Although the gas of inter- 
acting Bosons in ID is described by the exactly solvable 
Lieb-Liniger (LL) model the question of superfluid- 

ity is not readily answered from the exact solutions. The 
dimensionless interaction parameter 7 = g^m/{h^n) 
governs the crossover from the weakly-interacting quasi- 
condensate for 7 1 strongly-interacting Tonks- 

Girardeau gas [3j at 7 = 00, where is the coupling 
constant of the ID Bose gas ^ , n is the line density, and 
TO is the particle mass. The Landau criterion of superflu- 
idity [3 predicts a critical velocity Vc for the breakdown 
of dissipationless flow if no elementary excitations of mo- 
mentum p are accessible with energy below pvc to dissi- 
pate its energy. The excitation spectrum of the LL model 

contains umklapp excitations at finite momentum 2?™ 
and energies which tend to zero in a large system predict- 
ing a critical velocity of zero by the Landau criterion for 
any value of j [see the lower thin (blue) line in Fig. 
On the other hand, Luttinger-liquid theory and instanton 
techniques predict superfluidity in the LL model for 
sufficiently small 7. A resolution of this apparent para- 
dox lies in the probability of excitation by an infinitesimal 
external perturber which is given by the dynamic struc- 
ture factor (DSF) S{q,u!), the Fourier transform of the 
time-dependent density-density correlation function. A 
suppression of transitions in the vicinity of the umklapp 
excitation was found in Refs. 0, Il0l |. In this Letter, 
we report a calculation of the DSF, valid for large 7 as 
shown in Fig. We find that even at fairly large but 
finite 7 the low energy umklapp excitations are indeed 
suppressed while higher-energy excitations are enhanced. 
This is to be contrasted with the Tonks-Girardeau gas at 
7 = 00 where all excitations are equally well accessible 
and prevent superfiuidity. 

Experiments have recently made significant progress 
in confining Bose gases in ID probing the strongly- 
correlated Tonks-Girardeau regime by increasing interac- 
tions up to values of 7 « 5.5 and to effective values 
of 7off ~ 200 in an optical lattice 0| . In a recent exper- 
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FIG. 1: Excitation spectrum at 7 = 13. The upper and 
lower thin (blue) lines show the dispersions u!+{q) and uj-{q) 
of Eq. dJ, respectively, limiting the elementary excitations 
of the LL model. The dynamic structure factor S{q,uj) in the 
approximation Q is shown in shades of grey, between zero 
(white) and 1.0 N/{hqeY) (black). Here fcp = Tin and ep = 
h'^kl/{2m). The the presence of a 5-function contribution is 
indicated by red dots. 



iment jl4| . the zero momentum excitations of a ID Bose 
gas in an optical lattice have been measured by Bragg 
scattering, a technique that could also be used to mea- 
sure the DSF. 

Although the LL model is exactly solvable, it is no- 
toriously difficult to calculate the various correlation 
functions. Many results in limiting cases are summa- 
rized in the book [isj . but the full problem is not yet 
solved. Recently, progress has been made on various 
time- independent correlation functions 0, 0|. While 
a wealth of information is available for small 7 where 
Bogoliubov's perturbation theory can be applied, the 
strongly-interacting regime was hardly accessible as a 
systematic expansion in 7"^ was lacking. 

In this Letter we study the repulsively interacting ID 
Bose gas in a fermionic representation, introduced as fol- 
lows. The Bose-Fermi mapping of Girardeau 01 can be 
generalized to map the ID Bose gas with pair in- 
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teractions V{x) = gBS{x) to an equivalent interacting 
Fermi system with the same excitation energies and the 
same absolute values of the associated wave functions. 
Hence, any observable has the same value in both rep- 
resentations if it derives from a local operator in coor- 
dinate space, i.e. is a function of the density operator 
^(2^) = " ^i)- particular, this is the case 

for the DSF S{q,u!), which is the Fourier transform of 
the density-density correlation function , and ex- 

presses the probability to excite a particular excited state 
through a density perturbation 



S{q,uj) = 



'\pg\n)\H{noj-En + Eo), (1) 



where = X)i Gxp(— igx^) is the Fourier component of 
the density operator. In this mapping, large interactions 
between bosons correspond to small interactions between 
fermions, which allows for simple variational or perturba- 
tive treatment. However, the short-range interaction of 
Cheon and Shigehara |l8l | is highly singular for this pur- 
pose. It has been pointed out by Sen |20|| that the ground 
state energy functional can be derived variationally with 
the pair-interaction pseudopotential 



Vson(2;i,a:2) = -gFS"{xi - X2), 



(2) 



where gp = 2ti^ /{m^g-Q) is the coupling constant in the 
fermionic representation. This pseudopotential, however, 
is applicable only if the variational functions are contin- 
uous and vanish whenever two particle coordinates co- 
incide. This is the case for Slater determinants which 
we will use later to derive the Hartree-Fock (HF) and 
Random-Phase approximations (RFA) but not for the 
exact fermionic wavefunctions [l8j. Vson thus generates 
the first order correction to the ground-state energy easily 
whereas Sen found a rcnornializable divergence in second 
order perturbation theory. 

Recently, a more general pseudopotential was proposed 
by Girardeau and Olshanii (21|], which is also defined for 
discontinuous functions and generates the correct energy 
functional for arbitrary values of the coupling constant. 
We give the following useful representation of this pseu- 
dopotential in terms of an integral kernel: 
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Based on a pseudopotential Hamiltonian (0) or ((JJ, 
we can derive the HF approximation and the general- 
ized RPA giving access to the density-density correla- 
tion function and the static and dynamic structure fac- 
tors. The results should be valid at least to first order in 
7^^ and thus be useful near the Tonks-Girardcau limit 
of 7 = 00, where the interacting boson problem maps 
exactly to a free Fermi gas [J . The purpose of this work 
is thus twofold: Besides predicting the behavior of cor- 
relation fimctions of the ID Bose gas, we also test the 



validity and usefulness of the fermionic pseudopotential 
approach within the RPA. 

The HF approximation for the fermionic system is de- 
rived variationally in the standard way. We find that 
both pseudopotentials lO and ^ lead to exactly the 
same result: The HF equations for the single-particle or- 
bitals Lpj (x) take the usual form 



Fipj{x) = Sjipjix). 



(4) 



Unlike the case of Coulomb interactions where the Fock 
operator F is nonlocal with a local Hartree and a nonlocal 
exchange term wc find a purely local Fock operator 
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The first two terms on the right hand side are just the 
single-particle part of the Hamiltonian. The mean-field 
parts in the square bracket involve the single-particle 
density n{x) = ^-iid the derivative 

densities Vix) ~ ~iJ2j ^j'P*'j{^)'Pj{^) ^^'^ ^(•^) — 
n-j[<(5*(a;)(p"(a;) + 2ip*'^{x)ip'^{x)], reminiscent of mo- 
mentum and energy densities, respectively. Here, ip' ~ 
dip/dx and Uj is the fermionic occupation number Uj ~ 
l/{exp[(£j - fj.)/{kBT)] + 1}. At temperature T = 0, Uj 
defines the Fermi sea. 

For the homogeneous gas (T4xt = 0), the quantum 
number j is associated with the particle momentum, and 
the single-particle orbital solutions of Eq. arc plane 
waves with energy and effective mass 



(6) 



m* =m/ (1 - 2gFmn/h^) = m/{l - 47"^), (7) 

respectively. The HF approximation for the ground-state 
energy coincides with the first two terms of the large-7 
expansion of the exact ground-state energy in the LL 
model Q. 

The HF approximation permits us to calculate the lin- 
ear response function of time-dependent HF, also known 
as RPA with exchange or generalized RPA We 
calculate the DSF at zero temperature by the relation 
[l9| S{q,uj) = Im{x(9, — w — is)/''^} for a> > and 
S{q,uj) = for a; < from the dynamic polarizabil- 
ity xili'-^ + i^)- The function xili'-^ + determines 
the linear density response to an external field 0, 0| 
and can be obtained from the linearized equation of mo- 
tion of the density operator in the time-dependent HF 
approximation. In the Fourier representation of momen- 
tum and frequency, algebraic equations result which can 
be solved for x{q,to). In the thermodynamic limit, the 
sums over momenta are replaced by integrals. Wc find 
X(q,z) = X^°Kq,z)/{{l - 47-l)[i? - Dx^'^Kq^z)]}, with 
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z = oj + ie and + 0. For the DSF this yields for w > 



S{q,Lj) 



where 
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with the Tonks gas' Fermi wavenumber fcp = 
energy sf = fi^kp/{2m). The polarizabihty x*^"' 
1X2"' of the ideal ID Fermi gas with renormalized mass, 
is given by the real and imaginary parts 



X2 [q,^) = - 



Nm* 



■In 



Nnm* 



OJ 



I ±1, W_ < ±UJ < 

1 0, else, 



2hqkF 

where the dispersion relations 

LU±{q) = h\2kFq±q^\/{2m*) 



(12) 
(13) 

(14) 



border the continuum part of the accessible excitation 
spectrum made up from HF quasiparticle-quasihole exci- 
tations © , see Fig. ^ The (5- function part of the DSF 
© relating to discrete excitations of collective charac- 
ter lies outside of the continuum part and comes from 
possible zeros in the denominator of x{Qj^)- It is deter- 
mined by the solution ivo^q) of the transcendental equa- 
tion B — Dxi^ = 0. We have solved this equation in var- 
ious limits and found that at most one solution for ujq (q) 
may exist. The strength A(q) is given by the residue 
of X the pole z = uJo{q) — ie. For small 7"^, the 
strength A{q) ~ 2iV7exp(— 7g/fcF) is exponentially sup- 
pressed and possible solutions arc close to the dispersions 
cjg ~ ^±1 oc exp(— 7g/fcF). Due to this proximity of the 
discrete and continuous parts and expected smearing of 
discrete contributions beyond RPA, we may conjecture 
that the (5-function should be seen as part of the contin- 
uum, enhancing contributions near the border. In fact, 
we know from the exact solutions that energy spectrum 
is continuous 

At finite 7 we find a (5-function contribution below and 
close to LU- for small q whereas for large q there is a dis- 
crete contribution at energies larger than 10+ , see Fig. ^ 
In the limit q ^ 00 at finite 7, the 5 part completely 
determines the DSF as the continuum part vanishes, and 
asymptotically A ~ A^, and ujq ~ hq^ /{2m) becomes the 
free particle dispersion reminiscent of the DSF for the 
weakly interacting Bose gas in Bogoliubov theory Q. A 
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FIG. 2: DSF S{q, oj) as a function of a; at q = QAk^ [(a) and 
(b)] and q = 2fcF [(c) and (d)], and at interaction strength 7 = 
13.33 [(a) and (c)] and 7 = 40 [(b) and (d)]. The solid (black) 
line shows the RPA result ||HJ , the dashed (red) line shows the 
first order expansion in 7^^ 1151 1. and the dotted (blue) line 
shows the DSF of the Tonks limit (7 = 00) for comparison. 
While the first-order expansion always has a divergence to 
-l-oo near uj+, the RPA result shows an enhancement of low- 
energy excitation near U- at small q. Note the unphysical 
negative values of the DSF in the first order expansion near 
UJ-, in particular for the umklapp excitations in panel (c). 



more detailed analysis of the (5-function part will be given 
elsewhere |2^ . 

The branches lo± as shown in Fig. ^ ^^rs large 7 ap- 
proximations for the type I and II elementary excitation 
branches, respectively, introduced by Lieb In accor- 
dance with the exact results, both branches share the 
same slope at the origin and give rise to a single speed of 
sound given by Vs = dw± /dfc = UUf /rn* , which is the cor- 
rect first order expansion of Vs for large 7 . Note that 
the usual Bogoliubov perturbation theory for weakly in- 
teracting Bosons gives a similar expansion of Vg for small 
7 and the type I excitation branch. Type II excitations 
are not described in the Bogoliubov theory. 

Due to a logarithmic singularity in Xi'' the DSF van- 
ishes on the dispersion curves. For the Tonks gas at 
7 — > 00, the value of the DSF within these limits is inde- 
pendent of w and takes the value of N m / {2t:Ti^ qn) . The 
energy-dependence in the RPA for finite 7^^ is shown in 
Figs, n and |21 In particular, we see that the umklapp 
excitations aX q ~ 2/cf at small oj, prohibiting superfiuid- 
ity of the Tonks gas, are being suppressed. We find that 
S'(2/sF,(jj) in the RPA approaches zero as 1/In2(?iw/£F), 
in contrast to the results 0, ■, predicting a power-law 
dependence on oj for finite 7 based on a pseudoparticle- 
operator approach. 

The RPA result may be expanded in I/7, which yields 



S{q,uj)^ ^ fcp — 



87- 



ln/(g,t^) 
27 



0(7-^) (15) 



with f{q,uj) ~ |(ijj2 — (jj2)/(tj2^ — Li;2)|. This expansion is 
supposed to be consistent with first order perturbation 
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theory. However, it can assume negative values as seen 
m Fig. 01 although S{q, uj) is known to be strictly non- 
negative, a property that is fulfilled by our RPA result 
©• Close to the first order expansion has has a 
logarithmic singularity to +00 which may be a precursor 
of the dominance of Bogoliubov-like excitations in the 
DSF at small 7. In the RPA this effect is even more pro- 
nounced due to a strong and narrow peak of the DSF in 
the RPA near a;+ at large momenta as seen in Fig. |21 At 
finite gamma and for small momenta, however, the RPA 
predicts a peak near ui-, contrary to the first-order re- 
sult. Whether this effect is real or an artefact of the RPA 
is not obvious and may be decided by more accurate cal- 
culations or experiments. Spurious higher order terms in 
the RPA and an improved approximation scheme have 
been discussed in Ref. |2^. On the other hand. Roth 
and Burnett have recently observed a qualitatively sim- 
ilar effect in numerical calculations of the DSF of the 
Bose-Hubbard model ji^l • 

We have verified numerically that the /-sum rule mi = 
J ujS{q,Lu)du = Nn^q^/{2m) is fulfilled, which should 
be an exact statement for the RPA from a general the- 
orem The same follows also from the the large uj 
asymptotics of x{^.^^) — 2mi/{huj)'^ when x is analytic 
as a function of uj in the upper half complex plane 0. 
Our approximation breaks down when 7^8 for small 
values of momenta 5 ^ fcp, where imaginary poles of x 
appear. 

Our result (O defines approximations for the static 
structure factor S{q) = {h/N) J S{q,iLj)duj and the 
pair correlation function g{x) = 1 + J[S{q) — 
l\/(2'Kn) exTp(iqx)dx that will be discussed in detail else- 
where [231 . Here we only note that the approximation for 
S{q) is continuous and has the low-momentum expansion 
Siq) = q[l + Aj-^~qy{kh)]/{2kp) + 0{q^)h + 0ij-y 
We also find that g{x = 0) vanishes in first order of 7 ^ 
as consistent with Ref. , indicating once more the va- 
lidity of our results. 

Summarizing, we have derived variational approxima- 
tions that should yield valid expansions of various corre- 
lation functions of the ID Bose gas in I/7, making pre- 
dictions in a so-far unexplored regime. These approxima- 
tions prove consistent with known limits and sum rules 
and therefore establish the usefulness of the fermionic 
pseudopotentials ((SJ and ^ . We have also shown a rare 
example where a fully analytical calculation of the gen- 
eralized RPA can be carried out for a non-trivial many- 
body system. Our results have immediate relevance to 
current experimental endeavors to explore fermionization 
in the strongly-interacting ID Bose gas as they indicate 
qualitative deviation in the DSF from the Tonks gas limit 
already for very small 7^^. The more accessible region 
of intermediate values of 7, which would yield more in- 
sight into the superfluid properties, is beyond the scope 



of this study. Answers may come from higher order or 
density-functional-theory based approximations extend- 
ing the work of Ref. jl^, from numerical calculations, or 
from experimental measurements. Obvious extensions of 
the present study to finite temperature, inhomogeneous 
systems, and periodic lattice potentials are under way. 

The authors acknowledge discussions with Sandro 
Stringari who inspired their interest in this problem. 
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